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ABSTRACT 

We propose a method using the redshift dependence of the Alcock-Paczynski (AP) 
test and volume effect to measure the cosmic expansion history. The galaxy two-point 
correlation function as a function of angle, is measured at different redshifts. 

Assuming an incorrect cosmological model to convert galaxy redshifts to distances, the 
shape of ^(/i) appears anisotropic due to the AP effect, and the amplitude shifted by 
the change in comoving volume. Due to the redshift dependence of the AP and volume 
effect, both th e shape and amplitude of ^{fi) exhibit redshift dependence. Similar to 
iLi et al.l l|2014f) . we find the redshift-space distortions (RSD) caused by galaxy peculiar 
velocities, although significantly distorts exhibit much less redshift evolution 

compared to the AP and volume effects. By focusing on the redshift dependence of 
^(/i), we can correctly recover the cosmological parameters despite the contamination 
of RSD. The method is tested by using the Horizon Run 3 N-body simulation, from 
which we made a series of 1/8-sky mock surveys having 8 million physically self-bound 
halos and sampled to have roughly a uniform number density in z = 0 — 1.5. We find 
the AP effect results in tight, unbiased constraints on the density parameter and dark 
energy equation of state, with 68.3% CL intervals S^lm ^ 0.03 and 6w ~ 0.1, and the 
volume effect leads to much tighter constraints of S^lm ~ 0.007 and dw ^ 0.035. 


1 INTRODUCTION 

Since the discovery of c osmic acceleration dRiess et al.lll99^ : 
IPerlmutter et al]ll999l l. there have been many models pro¬ 
posed to explain it. These include modifications to general 
relativity (GR) and additional energy components to the 
constituents of the Universe. The simplest and currently 
favored solution is the additional field called dark energy 
which has negative pressure and is spatially homogeneous. 

To understand the nature of dark energy, we need to 
make precise measurements of the cosmic expansion history 
using cosmological observables. Two such observables are the 
angular diameter distance Da and the Hubble factor, H. If 
these can be measured at various redshifts, then tight con¬ 
straints can be placed on the parameters of the cosmological 
model, e.g., the matter density Ylm and the equation of state 
(EoS) of dark energy, w. 

Assuming an incorrect cosmological model for the co¬ 
ordinate transformation from redshift space to comoving 
space leaves residual geometric distortions. These distortion, 
known as the redshift-space distortions (RSD), are induced 
by the fact that distances along and perpendicular to the 
line of sight are fundamentally different. Measuring the ra¬ 


tio of galaxy clustering in the radial and transverse direc- 
tions, provides a probe of t he Alcock-Paczynski (AP) effect 
llAlcock fc Paczviiski|[l979l i. 


There have been several methods proposed that apply 
the AP test to the lar ge scale structure by measuring the 
clustering of galaxies (Ballingen IPeacock fc Heavens 199a : 


iMatsub ara fc Sutd 1 199a). s ymme tr y properties of galaxy 
pairs llMarinoiU-J^ Bu3 bOld : Ijmnings et_ajj ^ 


Bueno Belloso et al. I 20121'). and cosmic void s 1 RvderJ 


2011 


19951 : 


Lavaux fc Wandeltll2012l : ISutter et al.ll20lj '). Among them. 


the method of galaxy clustering has been widely used to 


constrain cosmological parameters dOutram et al 


Blake^f^jjjj20Uj;JChunngj^(Wana[j2012l: Reid^£3 


BeuHe£^t_^ 20ig;_Lindeiyetjjjj20jj: ISong et al.l 2014 : 


2004; 


2012 : 


Jeong et al.l 2014 : Lopez-Corredoira 20l4 ). 


The main caveat of this method is that, because the 
radial distances of galaxies are inferred from redshifts, AP 
tests are inevitably limited by RSD, which leads to ap- 
parent anisotropy even if the adopt ed c osmology is cor rect 
ilBallingeJIPeacock &; Heavens 1996l l. In iLi et ^ I (l2014h we 
proposed a new method utilizing the redshift dependence 
of AP effect to overcome the RSD problem, which uses the 
isotropy of the galaxy density gradient field. We found that 
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the redshift dependence of the anisotropy created by RSD is 
much less significant compared with the anisotropy caused 
by AP. Thus we focused on the redshift dependence of the 
gala:xy density gradient field, which is less affected by RSD, 
but still sensitive to cosmological parameters. 

The two-point correlation analysis is the most widely 
used method to study the large scale clust erings of galajc - 
ies. So, in this paper we revisit the topic of iLi et aP (l2014l l 
using the gala:xy two-point correlation function (2pCF). By 
investigating the redshift dependence of anisotropic galaxy 
clustering we can measure the AP effect despite of contami¬ 
nation from RSD. Moreover, when the redshift evolution of 
galaxy bias can be reliably modelled, we can measure the 
redshift evolution of volume effect from the amplitude of 
2pCF. The change of the comoving volume size is another 
consequence of a wrongly adopted cosmology, which has mo¬ 
tivated methods constraining cosmo logical parameters from 
number counting of g alaxy clusters llPress fc Shechtenl 19741 : 
IViana fc Liddlelll99^ and topology llPark fc Kimll2010l l. 

The outline of this paper proceeds as follows. In §2 we 
briefly review the nature and consequences of the AP effect 
and volume changes when performing coordinate transforms 
in a cosmological context. In §3 we describe the N-body 
simulations and mock galaxy catalogues that are used to test 
our methodology. In §4 and §5, we describe our new analysis 
method for quantifying the anisotropic clustering as well as 
proposing a way to deal with the RSD that are convolved 
with the AP distortion. Here we also present results of our 
optimised estimator. We conclude in §6. 


2 AP AND VOLUME EFFECT IN WRONGLY 
ASSUMED COSMOLOGIES 


In this section we briefly introduce the AP and volume effect 
in wrongly assumed cosmologies. A more detailed descrip¬ 
tion of the AP effect has been provided in lLi et ^ 

Suppose that we are probing the shape and volume of 
an object in the Universe. We measure its redshift span Az 
and angular size A9, then compute its sizes in the radial and 
transverse directions from the relations of 

Ary = Ar± = {1 + z)DA{z)Ae, (1) 

where H is the Hubble parameter. Da is the angular diam¬ 
eter distance. In the particular case of a flat universe with 
constant dark energy EoS, they take the forms of 


H{z) = Ho^Qma-3 + (1 - 


Da{z) 


1 , , 1 dz' 

1 + z^ ^ 1 + z Jo H{z') ’ 


( 2 ) 


where a = 1/(1 -|- z) is the cosmic scale factor, Hq is the 
present value of Hubble parameter and r{z} is the comoving 
distance. 

In case we adopted a wrong set of cosmological parame¬ 
ters in Equation (HU, the inferred Ary and Arx are wrong, 
resulting in distorted shape (AP effect) and wrongly esti¬ 
mated volume (volume effect). The degree of variations in 
shape and volume can be described by the following quan¬ 
tities 


[Ary/Arx] wrong _ [Da{z)H {z)]tiue 
[Ary/Arx] true [Da (2:)] wrong 


[Ar||(Ar_L)^] wrong Volumewrong [^a(^) /H(2)]\vrong /... 

[Ary(Arx)2]true Volumetrue [DA{zP/H(z)]tTue 

where “true” and “wrong” denote the values of quantities in 
the true cosmology and wrongly assumed cosmology. From 
the AP and volume effects, we can constrain Da{z)H{z) and 
Da{zY/H{z), respectively. 

The AP and volume effect due to wrongly assumed cos¬ 
mological parameters are shown in the upper panel of Fig¬ 
ure [T] Suppose that the true cosmology is a flat ACDM 
with present density parameter Q,m = 0.26 and standard 
dark energy EoS w = — 1. If we were to distribute four per¬ 
fect squares at various distances from 500 Mpc/h to 3,000 
Mpc/h, and an observer located at the origin were to mea¬ 
sure their redshifts and compute their positions and shapes 
using redshift-distance relations of two incorrect cosmolo¬ 
gies: 

(i) = 0.41, w =-1.3, 

(ii) = 0.11, w = -0.7, 

the shapes of the squares appear distorted (AP effect), and 
their volumes are changed (volume effect). In the cosmologi¬ 
cal model (ii) with Dm, = 0.11, w = —0.7, we see a stretch of 
the shape in the line of sight (LOS) direction (hereafter “LOS 
shape stretch”) and magnification of the volume (hereafter 
“volume magnification”), while in the model with Dm, = 0.41, 
w = —1.3, we see opposite effects of LOS shape compression 
and volume shrinkage. 

In the lower panel of Figure [T] we plot the degree of 
variations in shape and volume as functions of redshift. In 
cosmology (i), both quantities have values less than 1, indi¬ 
cating LOS shape compression and volume shrink. The effect 
in cosmology (ii) is slightly more subtle. At low redshift, the 
effect of dark energy is important, and there is LOS shape 
compression and volume reduction due to the quintessence 
like dark energy EoS. However, at higher redshift the role 
of dark matter is more important, and we see LOS shape 
stretch and volume magnification due to the small Q.m- 

More importantly. Figure [T] highlights the redshift de¬ 
pendence of the AP and volume effects. For example, in the 
cosmology with = 0.41, w = —1.3, both the LOS shape 
stretch and volume magnification become more significant 
with increasing redshift. In the cosmology with Dm, = 0.11, 
w = —0.7, not only do the magnitudes of the effects evolve 
with redshift, but there is also a turnover from LOS shape 
compression and volume shrink at lower redshift to LOS 
shape stretch and volume magnification at higher redshift. 


3 MOGK DATA 

We test the methodology using mock surveys constructed 
from one of the Horizon Run simulations, HR3. The Hori¬ 
zon Run simulations are a suite of large volume N-body 
simulations that have resolutions and volumes capable of 
reproducing the observational statistics of many current 
major redshift surveys like the Sloan Digital Sky Survey’s 

i SDSS) Baryon Oscillation Spectroscopi c Survey (BOSS) etc 
Park et alK2005l : lKim et al.ll2009l . 1201 il l. HRS adopts a flat- 
space ACDM cosmology with the WMAP 5 year parameters 
Dm = 0.26, Hq = 72 km/s/Mpc, = 0.96 and as = 0.79 
l|Komatsu et al.ll201lll . The simulation was made in a cube 
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Figure 1. The redshift dependence of AP and volume effect in two wrongly assumed cosmologies flrn = 0.41, w = —1.3 and Qrn = 0.11, 
w = —0.7, assuming a true cosmology of Q,rn = 0.26, w = —1. Upper panel shows the apparent distortion of four perfect squares, 
measured by an observer located at the origin. The apparently distorted shapes are plotted in red solid lines. The underlying true shapes 
are plotted in blue dashed lines. Lower panel shows the evolution of Equations 0 and 0. In our mock surveys we split the samples at 
^ = 0.3, 0.6, 0.9 and 1.2, as marked by the vertical lines. 


of volume (10.815 h ^Gpc)® using 7120® particles with par¬ 
ticle mass of 1.25 x 10 ^®/i“®Mq . 

The simulations were integrated from z = 27 and 
reached ^ = 0 after making N^tep = 600 timesteps. Dark 
matter halos were identified using the Friend-of-Friend al¬ 
gorithm with the linking length of 0.2 times the mean par¬ 
ticle separation. Then the physically self-bound (PSB) sub¬ 
halos that are gravitationally sel f-bound and tidally stable 
are identified dKim &; Park|[200^ ). The PSB halo finder is a 
group finding algorithm which can efficiently identify halos 
located even in crowded regions. This method combines two 
physical criteria such as the tidal radius of a halo and the 
total energy of each particle to find member particles. 

An all-sky, very deep light cone survey reaching redshift 
z = 4.3 was made by placing an observer located at the 
center of the box. The comoving positions and velocities of 
all CDM particles are saved as they cross the past light cone 
and PSB subhalos are identified from this particle data. 

To match the observations of rec ent LRG surveys 
llGott et al.ll200i . l2009l : IChoi et al.ll2010l) . a volume-limited 
sample of halos with constant number density of 3 x 
10“^(/i“®Mpc)“® is selected by imposing a minimum 
halo mass limit varying along with redshift, from 6 x 


at high redshift to 1.5 x at low 

redshift. To do this, the whole sample is split into many 
thin shells whose radial spans are 30 Mpc/h. In each shell, 
a minimal mass cut is set such that the number density is 
h = 3 X 10“^(/i“®Mpc)“®. These thin shells are then merged 
together to build up the whole sample. Note that all these 
steps are done before converting the sample to other cos¬ 
mologies, so the constant number density of halos is only 
true in the correct cosmology. The group velocity of each 
PSB dark matter halo (or our mock “galaxy”) is used to 
perturb the redshift of the halo by 


Az = (1 + z) 


ulos 

c 


(5) 


where ulos is the LOS component of the halo proper veloc¬ 
ity. 

We divide the all-sky survey sample into eight equal sky 
area subsamples and impose the redshift range z = 0 — 1.5. 
This mock data will be re levant for future galaxy spectro¬ 
scopic surveys (e.g. DESI iLevi et ^ liH). 
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4 METHODOLOGY 


We probe the effects discussed in ^ using the 2pCF. The 
2pCF is a mature statistic in cosmology and its optimal 
estimation considers minimal variance while dealing with 
complicated masks and selection functions. The most com¬ 
monly a dopted formulation is that of the Landy-Szalay es¬ 
timator jLandv fc Szalavl|l993l l. 




DD - 2DR + RR 
RR 


( 6 ) 


where DD is the number of galaxy-galaxy pairs, DR the 
number of galaxy-random pairs, and RR is the number of 
random-random pairs, all separated by a distance defined 
by s ± As and p ± Ap, where s is the distance between 
the pair and p = cos{6), with 0 being the angle between 
the line joining the pair of galaxies and the LOS direction. 
This statistic therefore captures the radial anisotropy of the 
clustering signal. 

The random point catalogue constitutes an unclustered 
but observationally representative sample of our mock sur¬ 
veys. To reduce the statistical variance of the estimator we 
use 15 times as many randoms as we have galaxies. 

To probe the anisotropy, ^ is measured at different di¬ 
rections, and integrated over the interval As = Smax — Smin. 
We evaluate 


/ ■Smax 

Gs,n)ds. (7) 

min 

We limit the integral at both small and large scales. At small 
scales the shape of (.As ia) is seriously distorted by the fin¬ 
ger of god (FoG) effect ll.lackson Ill972l l. and the distortion is 
more significant at lower redshift where structure undergoes 
more non-linear growth. This introduces redshift evolution 
in ^As(m) which is rather difficult to model. As will be men¬ 
tioned in §5, since FoG is particularly strong near the LOS 
direction, in our analysis we also impose a cut /r > 0.1 to 
further suppress its contamination. At large scales the mea¬ 
surement is dominated by noise due to poor statistics. We 
find Smin = 6—10 Mpc/h and Smax = 40 — 70 Mpc/h are 
reasonable regions which can provide consistent results of 
2pCF and tight, unbiased constraints on cosmological pa¬ 
rameters. In our analysis we choose Smin = 6 Mpc/h and 
Smax ■— 50 Mpc/h. 

As an example Figure [2] shows how the 2pCF is affected 
by AP and volume effects in the Q„i = 0.41, w = —1.3 
cosmology. In choosing incorrect cosmological parameters, 
we expect the 2pGF to be influenced in three ways. First, 
as a result of the AP effect, structures appear compressed 
in the radial direction. This induces a nonuniform variation 
in ^As{n) ^ a- function of angle. Second, as a result of the 
volume effect, the sizes of structures shrink. For example, 
a structure whose original size is so = 50 Mpc/h will show 
up with a size si < 50 Mpc/h. As a result, the amplitude 
of ^Asin) changes. Finally, as another consequence of the 
volume effect, in the wrong cosmology structures on larger 
scales enter the statistics. For example, halos within the blue 
solid box are not considered in the correct cosmology, but 
they are taken into consideration in the wrong cosmology, as 
shown by the black dashed box. This also results in a change 
in the amplitude of ^As(/r) since the binning in s,/i-space 
will be inconsistent between different cosmological models. 
The combined effects of choosing an incorrect cosmology on 


shear and volume have been noted bv IPark fc KimI ll201(]ll . 
who used the volume effects measured by the genus statistic 
to constrain the expansion history of the Universe. 


5 RESULT 

In this section we present the measurements of 2pGF from 
our mock data. Figure [3] shows the ^as measured from 
HRS mock surveys, adopting the correct cosmology (left), 
the Qm = 0.11, w = —0.7 cosmology (middle), and the 
Um = 0.41, w = —1.3 cosmology (right). To study the red¬ 
shift evolution, we divide the redshift range 2 = 0— 1.5 
into five equal-width redshift binsQ. To measure anisotropy, 
we further divide the full angular range /r = 0 — 1 into 10 
equal-width bins. So we have 

^Aa(zi, Mj) = Cas in the i — th redshift bin, j — th p bin(8) 

where i = 1, 2,..., 5, j = 1,2, ...10. Measurements without 
and with considering RSD effect are plotted in solid and 
dashed lines, respectively. 

The result for the correct cosmology is plotted in the 
upper left panel of Figure 0 In the absence of RSD, we ob¬ 
tain flat curves of ^Aa{n) in nil redshift bins, with the am¬ 
plitude slightly different from one redshift bin to another. 
This difference can arise from two sources; (a) The growth 
of clustering with the decreasing of redshift. (b) The redshift 
evolution of the bias of halos having the same comoving den¬ 
sity. The result is significantly changed when we include the 
RSD effect. Near the LOS direction (p ^ 1), structures are 
compressed due to the Kaiser effect, so the value of ^Aa{D 
is smaller compared with measurements near the tangential 
directior0. But it should be noted that the shape of ^Aain) 
is nearly the same at all redshifts, indicating the small red¬ 
shift dependence of the RSD effect. It is this observation that 
makes our method both feasible and statistically powerful. 
Even though the 2pGF becomes very anisotropic in redshift 
space, the anisotropy due to RSD does not change much as 
a function of redshift and its redshift-dependence is dom¬ 
inated by the geometric effects introduced by the adopted 
cosmology. 

The results of the Qm = 0.11, w = —0.7 and Qm = 0.41, 
w = —1.3 cosmologies are plotted in the middle and right 
panels, respectively. We can see that ^as is significantly al¬ 
tered by the volume and AP effects. In the flm = 0.41, 
w = —1.3 cosmology, the shrinkage of comoving volume sup¬ 
presses the amplitude of ^as(m), and the LOS shape com¬ 
pression of structures results in a suppression of amplitude 
in the LOS direction compared with the tangential. Both 
effects become increasingly more significant at higher red¬ 
shift. In the flm = 0.11, w = —0.7 cosmology, the comoving 
volume shrinks at 2 < 0.6 and enlarges at 2 > 0.6. Thus, 
compared to the result of correct cosmology (top left panel), 
in this cosmology the 2pCF measured at 2 = 0.3 — 0.6 has 
slightly suppressed amplitude, while the other three curves 
measured at 2 > 0.6 have enhanced amplitudes. Also, due 


^ We do not show the result of the first redshift bin, which is 
noisy due to poor statistics. 

^ The FoG effect will enhance ^Aa{G in the LOS direction. It 
does not significantly show up in our figures since we impose the 
cut Smin = 6 Mpc/h. 
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Figure 2. The influence of the AP and volume effect on the 2pCF statistics, in the wrongly assumed cosmology Qrn = 0.41, w = —1.3. 
Scatter points are a 200 Mpc/h x 200 Mpc/h slice of halos drawn from HRS, with thickness 10 Mpc/h and mean redshift 2 : = 0.43. Their 
positions are measured by an observer located at the origin, with the LOS direction marked by the green arrow. Blue empty circles are 
their true positions in the correct cosmology, while red filled circles are the apparent positions in the wrong cosmology. Blue solid circles 
are the “2-point CF rings” (i.e. the geometrical bins used to evaluate the clustering) with radii 6 and 50 Mpc/h, centered on a galaxy 
with location x,y=1200 Mpc/h. Red solid circles are the “exactly mapped” rings (bins) in the incorrect cosmology, with its size shrunk 
and shape distorted. Black dashed circles are the actually adopted rings when we measure the 2pCF in the wrong cosmology. We expect 
the result of 2pCF statistics to be affected in three ways. First, the shape of structures should appear compressed in the radial direction, 
introducing anisotropy in Secondly, the sizes of structures should appear shrunk due to the volume effect, which shifts the amplitude 
of ^ . Finally, as another consequence of volume effect, structures on larger scales enter the statistics. For example, halos within blue 
solid box are not considered in the correct cosmology, but will be considered in the wrong cosmology, as shown by the black dashed box. 
This also shifts the amplitude of 


to the LOS shape stretch at 2 > 0.3, in all curves the ampli¬ 
tude is relatively enhanced in the LOS direction. All these 
effects are more significant at earlier times. 

In real observational data the redshift evolution of the 
bias of observed galaxies is difficult to model. Thus to miti¬ 
gate this systematic uncertainty we wish to rely on the shape 
of rather than its amplitude. In the lower panel of 

Figure [3l we show the normalized ^As(fi) in each redshift 
bin, defined as 


iAsin) 




(9) 


It should be noted that the normalization in Equation ((Oj 
is made by the correlation function in each corresponding 
redshift bin, not by the global quantity. This procedure is 
critically important since it effectively suppresses the cos¬ 
mic variance in the correlation function from one redshift 
bin to another. When the RSD effect is considered, in the 
correct cosmology ^as at different redshifts are almost iden¬ 


tical to each other, while in the wrong cosmologies we see 
clear redshift evolution. 

Overall, the effect of RSD on the 2pCF is large but its 
redshift dependence is small. The correct cosmology corre¬ 
sponds to the case with the lowest change of ^as with 2 . 
Even with RSD, we can still correctly determine the true 
cosmology by using the relative change of ^as with redshift. 
Based on this fact, we define our ^is follows 

[iAsjzi, flj) - ^Asjzs, 

i=l j = l 

The 2pCF measured in the 1-4 redshift bins are compared 
(or normalized) to the measurement in the last redshift bin. 
This will prefer minimal shape change over the redshift 
range, with little of no weight given to the amplitude of the 
clustering statistic. 

Figure [ 3 ] shows that the redshift evolution of the RSD 
effect, although small, still results in visible redshift evolu¬ 
tion in the amplitudes and shapes of the curves. Also, the 
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Figure 3. The 2pCF measured in four redshift bins, in the correct cosmology (left) and two wrongly assumed cosmologies (middle: 

= 0.11, w = —0.7; right: Qrn = 0.41, w = —1.3). The clustering signal is measured as a function of 1 — /r, where fi = cos{6) and 9 
is the angle between the LOS and the vector joining the pair of galaxies. Dashed and solid lines show the results with and without the 
RSD effect, respectively. Upper panel: In the wrongly assumed cosmologies, we observe a clear change in the amplitudes and shapes of 5 
due to the volume and AP effect. Additionally, due to the redshift dependence of volume and AP effect, the amplitudes and shapes in 
the four redshift bins are different from each other. Lower panel: The same as the upper panel, except that the amplitudes of curves are 
normalized to 1. 


amplitude of the 2pCF may evolve due to the growth of 
clustering and the r edshift evolution of selection effect. So, 
similar to iLi et ahl (|2014h . we further correct the residual 
RSD effect, i.e., the following quantity is computed in the 
correct cosmology and subtracted from our results, 

A^As = ^As.With RSD — Cas.No RSD- (11) 

Although this correction factor may have cosmological de¬ 
pendence, we assume that this will not introduce significant 
systematic to our methodology. To avoid the contamination 
from the FoG effect we further impose a cut of /r > 0.1 


(i.e. j > 1). We also take into consideration the correlation 
between measured in different directions. So the 

function becomes 

4 10 10 

PGi: fljL, ( 12 ) 

i=l ji=2j2=2 

where p{zi,^j) takes the form 

p(Zi, fj,j) = ^As{Zi, fj^j) — ^As(-2^5) f^j) ~ 

fij) — A^As(' 2;5 7 Mj)]) 


(13) 
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and Covi is the covariance matrix of p{zi,fij)- In testing our 
methodology, we divide the full sky sample into 48 equal sky 
coverage subsamples and construct a measurement covari¬ 
ance matrix. 

Figure |4] shows the likelihood contours obtained from 
our mock survey, which is a 1 / 8 -sky survey with constant 
number density h = 3 x 10“"^(/i~^Mpc)“® in the redshift 
range 0-1.5. The gray contours on the left show the result 
when we normalize the amplitude of ^Asiu) and just utilize 
the information encoded in the shape (i.e., we focus only 
on the AP information and ignore the volume effect). We 
obtain tight and unbiased constraint on flm and w, with 
68.3% CL intervals Silrn 0.03 and Sw ~ 0.1. We find that 
tlm and w are positively degenerated with each other, which 
is expectable. For instance, reducing Qrn and having a more 
phantom-like dark energy produce similar influences on the 
expansion history of the Univ erse. Compared with the result 
of our gradient field method |Li et al.l l|2014h . the shapes of 
contours are very similar, and constraint is a little tighter. 

The right hand plot of Figure |4] shows the result when 
we utilize the volume effect. To do that, in Equation (|13l) 
we adopt the non-normalized ^as instead of the normalized 
one ^As 0. In that case, constraints become much tighter, 
with S^lrn ~ 0.007 and Sw ~ 0.035. Also, the direction of 
degeneracy changes and is very different from mainstream 
techniques of CMB, SNIa and BAO, meaning that combin¬ 
ing our method with these techniques can significantly im¬ 
prove the constraint. To implement it in real observational 
cases, it is necessary to model the evolution of the clustering 
amplitude for the observed galaxies. 


6 CONCLUSION 

We have presented a new anisotropic clustering statistic that 
can probe the cosmic expansion history. We measure the in¬ 
tegrated 2pCF, ^As{u) = as a function of 

direction /r. The amplitude of ^as(m) is affected by the vol¬ 
ume effect, and the shape is affected by the AP effect. Due 
to the redshift dependence of the volume and AP effects, in 
wrongly adopted cosmologies there are redshift evolutions of 
the amplitude and shape. The RSD effect due to galaxy pe¬ 
culiar velocities, although having a strong effect on ^As{n)j 
does not exhibit significant redshift evolution. Thus by fo¬ 
cusing on the redshift dependence of ^as (u) > we are able to 
derive accurate and unbiased estimates of cosmological pa¬ 
rameters in spite of contamination induced by RSD. We are 
only assuming that the RSD modelling is accurately cap¬ 
tured by the simulations. _ _ 

The concept of this paper is similar to iLi et ^ 
where the redshift dependence of the AP effect is mea¬ 
sured from the anisotropy in the galaxy density gradient 


^ Different from the normalized case, in the non-normalized 
case not only do we correct the RSD effect, but we also 
correct the redshift evolution of the amplitude of 2pCF. As 
shown by the top left pannel of Figure [3] the amplitude of 
2pCF can evolve with redshift even in the no RSD case, 
mainly due to the growth of structure and the selection bias. 
So in the non-normalized case Equation mj shall be mod¬ 
ified as A^As{zi) = [?As,With RSaC^i) — ?As,No RSD(^i)] + 
[Sas.No RSoC^^i) — Cas.No RSd(^i)] • 


field. However, in this paper we choose a different statistical 
method, i.e. the 2pCF. They differ from each other in several 
aspects. 1) Using the 2pCF method it is more convenient to 
choose the scales we investigate. 2) The advantage of the 
density gradient field method is that, it allows us to utilize 
the information on small scales of -^10 Mpc/h (depending 
on the scale of smoothing). 3) In the 2pCF method we are 
able to probe the volume effect, which is not possible for the 
galaxy density field method. 4) The 2pCF is a mature statis¬ 
tic in cosmology and its optimal estimation and statistical 
properties are well understood. 

The volume effect, which causes redshift evolution in the 
amplitude of 2pCF, leads to very tight constraint on cosmo¬ 
logical parameters. But it suffers from systematic effects of 
growth of clustering and the variation of galaxy sample with 
redshift. It would be great if one can reliably model these 
two effects and utilize the volume effect. In case that the 
systematic effect can not be correctly modelled, one can fo¬ 
cus on the AP effect by normalizing the amplitude of ^as(m) 
and just investigating the redshift evolution of the shape. 

In our analysis we consider the case of a fiat uni¬ 
verse with constant dark energy EoS w. However, in general 
one can choose an arbitrary cosmological model and repeat 
our analysis to constrain the model parameters. Further¬ 
more, since our method focuses on the redshift evolution of 
the 2pCF, it may have advantages in constraining models 
in which the cosmic expansion history displays dynamical 
features. For the popular Chevallier- Polarski-Linder (CPL) 
param etrization w — wo + Wajf^ l|Chevallier fc Polarskil 
I2OOII : [LindeillioO^ L the dark energy EoS evolves monotoni- 
cally from uio at 2 = 0 to wq + Wa at high redshift. Using a 
LSS sample covering a wide redshift range should be suffi¬ 
cient to constrain the model parameters through the effect 



a transition of w could happen in a very short period, which 
is difficult to be detected by common techniques. For our 
method, we can reduce the redshift bin size down the a size 
similar to the typical width of a transition epoch to probe 
the detailed feature of cosmic expansion history and effec¬ 
tively constrain the model parameters. 

The feasibility of our method crucially depends on the 
fact that the redshift dependence of the RSD effect is small. 
One may suspect that this is a special consequence of the 
constant number density sample. While an exhaustive inves¬ 
tigation of this issue is too complicate for this short paper, 
we perform a simple check to show that the weak redshift 
dependence of the RSD effects is fairly universal. From the 
HR3 lightcone data, we generated a series of samples with 
a constant mass cut, and measured the root mean square 
(RMS) displacement of halo comoving distance as a function 
of redshift. The degree of the RSD effects shall be statisti¬ 
cally directly proportional to the RMS displacement. The 
result is plotted in Figure [S] We find that the total varia¬ 
tion of the RMS displacement between 2 = 0 and 1.5 is only 
about 13%, which can explain why the redshift dependence 
of the RSD effects is weak. We also find that, within the wide 
mass cut range from 5 x to 2 x lO^®/i“^M 0 , the 

difference in the RMS displacement is < 1%. This suggests 
that the redshift evolution of RSD effect shall be fairly in¬ 
sensitive to the halo mass cut and number density of the 
sample. 




















8 Xiao-Dong Li, Changhom Park, Cristiano G. Sabiu and Juhan Kim 



LI 


m 



Figure 4. Left: Expected likelihood contours (68.3%, 95.4%, and 99.7%) in the Qrn — w plane, obtained from a 1/8-sky, z < 1.5 mock 
survey with a constant galaxy number density of n = 3 X 10“^(h“^Mpc)“®. We achieve unbiased constraints with 5w ~ 0.1 and 
SQrn ~ 0.03 by comparing the shapes of 5 as(m) measured in different redshift bins. Right: Here we use the unnormalized 5 as(m)i which 
is sensitive to the volume change and thus provides much tighter constraints. Although to use this in practice would mean overcoming 
some observational systematic uncertainties like galaxy evolution and selection bias. 



Figure 5. The RMS displacement of halo comoving distance as a function of redshift, for samples with different cuts of halo mass. Over 
the redshift range of interest for our study, the peculiar velocity is varying smoothly with a maximal deviation of 13%. Within the wide 
mass cut range from 5 X lO^^/i“^M 0 to 2 X lO^®h“^M 0 , the difference in the RMS displacement is < 1%, suggesting that the redshift 
evolution of RSD effect shall be fairly insensitive to the halo mass cut and number density of the sample. 


When dealing with real observational data, it will be im¬ 
portant to accurately model the gala:xy clustering to remove 
the residual RSD effects on the 2pCF. It will also require 
the handling of various observation effects such as survey 
geometry, fiber collisions, etc. We will report the results of 
such investigations in forthcoming studies. 
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